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ABSTRACT 

Observational evidence has been mounting over the past decade that at least some luminous (~ 2L*) galaxies 
have formed nearly all of their stars within a short period of time only 1-2 x 10 9 years after the Big Bang. 
These are examples of the first major episodes of star formation in the Universe and provide insights into the 
formation of the earliest massive galaxies. We have examined in detail the stellar populations of six z ~ 1.5 
galaxies that appear to be passively evolving, using both ground and space-based photometry covering rest- 
frame UV to visible wavelengths. In addition, we have obtained medium-resolution spectroscopy for five of 
the six galaxies, covering the rest-frame UV portion of the spectrum. Spectral synthesis modeling for four of 
these galaxies favors a single burst of star formation more than 1 Gyr before the observed epoch. The other two 
exhibit slightly younger ages with a higher dust content and evidence for a small contribution from either recent 
star formation or active nuclei. The implied formation redshifts for the oldest of these sources are consistent 
with previous studies of passive galaxies at high redshift, and improved stellar modeling has shown these results 
to be quite robust. It now seems clear that any valid galaxy formation scenario must be able to account for these 
massive (~ 2 x 10 11 M Q ) galaxies at very early times in the Universe. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high redshift 



1. INTRODUCTION 

Understanding the formation histories of the most mas- 
sive galaxies has remained an important topic shaping both 
observational and theoretical studies of galaxy formation 
and evolution. The basic structure of a cold dark matter 
(CDM) cosmology requires a bottom-up formation mecha- 
nism with a late-time build-up of the most massive systems. 
Early semi-analytic models (SAMs) based on this paradigm 
predicted that massive galaxies did not assemble the ma- 
jority of their present-day mass until 1 < z < 2 (e.g ., 
iThomas & Kauffmannl [19991: iKauffmann & Haehnelti 2000). 
Evid ence for cosm ic "downsizing" ("e.g.. ICowie et al. 1996; 
lJuneau et al.ll2005h however, appears to be inconsistent with 
the results generated by these early hierarchical models. 
More recently, the additi on of feedback from active galactic 



nuclei (AGN) into SAMs dGranato et al. 



2006; [Bower et al.l 120061: iDe Lucia et ail 



2004 iCrotonetail 



2006) has provided 



a means to shut off star formation at earlier times and allow 
galaxies to evolve onto the red sequence much more quickly, 

1 Some of the data presented herein were obtained at the W. M. Keck Ob- 
servatory, which is operated as a scientific partnership among the California 
Institute of Technology, the University of California, and the National Aero- 
nautics and Space Administration. The observatory was made possible by the 
generous financial support of the W. M. Keck Foundation. Results are also 
based in part on data collected at Subaru Telescope, which is operated by 
the National Astronomical Observatory of Japan, and on observations made 
with the NASA/ESA Hubble Space Telescope, obtained at the Space Tele- 
scope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program # GO-10418. 

Visiting Astronomer at the Infrared Telescope Facility, which is operated 
by the University of Hawaii under Cooperative Agreement no. NCC 5-538 
with the National Aeronautics and Space Administration, Science Mission 
Directorate, Planetary Astronomy Program. 

3 Also at Cerro Tololo Interamerican Observatory, Casilla 603, La Serena, 
Chile. 



thereby producing massive "red and dead" galaxies at higher 
redshifts. The underlying hierarchical growth of structure, 
however, remains a fundamental assum ption throughout. I n 
the AGN feedback scenario proposed by lCroton et al.1 (120061) . 
the "quasar mode" operates between 2 < z < 4 during which 
time the central black hole is built up, until the more quiescent 
"radio mode" takes over, quenching star formation by a slow 
heating of the surrounding halo gas. 

Observational constraints have largely driven the need for 
these corrections to SAMs and remain a key aspect for contin- 
ued understanding of the processes of galaxy formation. One 
method for testing predictions from these models is to identify 
galaxies containing the oldest stellar populations at high red- 
shift. These serve as examples of the first major episodes of 
star formation in the universe and provide laboratories where 
we can study these processes in an essentially pure and un- 
confused environment. 

Several recent studies have discovered "old galaxies" (OGs; 
by which we refer to the age of the dominant stellar pop- 
ulation and not necessarily the age since galaxy assembly) 
among extremely red objects (EROs) defined to have R—K > 
5 - 6 dLiu et al.ll2000t iDaddi et alJi2000. 2005; Cimatti et al.1 



20021 12004 live et alj|2003t lYan et al J 12004 iFu et al J 120051 



iyj 

Stockton et al. 2004, 2006, 2007; Stockton & McGrath 2007). 



In order to discriminate between old stellar populations and 
dusty starburst galax ies which contribute roughly equal ly to 
the ERO population (iCimatti et al.ll2002t lYan et aLll2004l) . the 
strength of the 4000 A break can be used as a simple diagnos- 
tic feature. Old stellar populations (> 1 Gyr) show a sharp in- 
flection at 4000 A due to metal absorption in the atmospheres 
of late-type stars, whereas dust-reddened starbursts show a 
spectral slope that is much less inclined near 4000 A. Using 
the strength of this break, an approximate age can be deter- 
mined for the galaxy, which can be used to extrapolate back to 
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FIG. 1 . — Finding charts for the six OGs discussed in this paper. For each field, the quasar is centered in the image. North is up and East is to the left. 



the formation redshift. The results of these studies imply that 
a population of massive galaxies existed at very early times 
in the universe (zf > 3) and that star formation must have 
proceeded extremely rapidly in these early massive galaxies. 
These conclusions are consistent with "fossil" studies of lo- 
cal massive elliptical galaxies, which conclude that the stellar 
populations in s uch galaxies formed early and over a short 
time period (e.g., Thomas et al. 2005; Nelan et al. 2005|). 

We present a sample of six z ~ 1.5 OGs that have both 
R-K 1 > 6, and J-K' ~ 2. Finding charts for these OGs 
are given in Fig. [T] We have obtained broad-band photom- 
etry covering rest-frame UV to visible wavelengths, as well 
as medium resolution spectroscopy covering the rest-frame 
UV portion of the spectrum for five of the six sources. The 
spectra provide a more precise measure of the galaxy redshifts 
than those determined from the photometry alone, and these 
redshifts, in turn, allow a more accurate determination of the 
age of the dominant stellar populations. The spectra are also 
more sensitive to age diagnostics, such as the Mg II doublet at 
2800 A, as well as to star formation indicators, such as [O II], 
which are indistinguishable from the underlying continuum 
in the broad-band photometry. We examine the stellar popu- 
lations of these OGs by comparing our observations with the 
most recent population synthesis models, and place them in 
context with current semi-analytic models of galaxy forma- 
tion. Throughout this paper we assume a flat cosmology with 
H = 71 km s- 1 Mpc-\ Q M = 0.27 and Q A = 0.73. 



2. OBSERVATIONS 



Our sample of six OGs at z ~ 1.5 were found in the fields 
of radio-loud quasars. Radio quasar fields were selected based 
on both the knowledge that radio sources tend to pinpoint 
more overdense regions of sp ace than their optically selected 
radio-quiet counterparts (e. g., Best 200 0t ICroom et al!| | 2001b 
iBest et alj|200l iBarr et alj[2003l; iRottgerin g et al.ll2005l) and 
the assumption that star formation would proceed most 
rapidly in the most overdense regions of space, therefore en- 
abling the detection of evolved galaxies at high redshift. De- 
tails of the pho tom etric selection criteria can be found in 
IStocktonl (1200 ll) and lStockton et al.l (l2006h : however, we will 
provide a brief summary here. 

Fields were imaged in the if '-band, and we looked for ob- 
jects with magnitudes in the range 17 < K' < 19 expected 
for a galaxy at z ~ 1.5 that would passively evolve into a 
~ l-5£* galaxy by the present day. J-band photometry was 
then obtained to search for all objects with J — K' ~ 2, ex- 
pected for old stellar populations (> 1 Gyr) in this redshift 
range. Finally i?-band imaging was acquired to constrain the 
blue ward sid e of the 4000 A-break and look for objects with 
R-K' > 6. IStockton et alj (120061) show that these combined 
R—K' > 6 and J —K' ~ 2 criteria are sufficient to discrimi- 
nate between heavily dust-reddened star-forming galaxies and 
passively evolving old stellar populations at z ~ 1.5 (cf., their 
Figure 1). Once these criteria were met, we then obtained 
additional imaging at blue and intermediate wavelengths and 
deeper imaging in the near-IR in order to reduce photometric 
errors and better discriminate between model stellar popula- 
tions. 

This selection method differs in two important ways from 
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that of iDaddi etall d2004l) who use the BzK criterion to 
find massive galaxies in the Hubble deep fields. First, selec- 
tion around radio sources allows us to select specific redshift 
ranges where the diagnostics from broad-band filters gives the 
cleanest separation between dusty starbursts and old stellar 
populations. Specifically, for z ~ 1.5, the 4000 A-break falls 
between / and J filters. Second, by selecting sources in the 
infrared and using the ERO R — K' > 6 criterion (as opposed 
to requiring a S-band detection) we are able to reach the faint 
limiting magnitudes required to avoid bias against galaxies 
with little or no blue flux, thereby enabling us to pick out the 
oldest stellar populations more efficiently. This method also 
differs from the distant red galaxy (DRG) selection criterio n 
of J - K > 2.3 dFranx et al.ll2003t Ivan Dokkum et al.ll2003l) . 
which is primarily used to select galaxies at z > 2 without 
regard for the age of the dominant stellar population. 

2.1. Optical and Near-IR Ground-based Imaging 

Near-IR observations at J and K 1 were obtained on the 
NASA Infrared T elescope Facility (IRTF) with NSFCam 
dShure et al.lll994l) for all of the fields. Deeper follow-up J 
and K 1 photometry w ere obtained on Subaru with CISCO 
dMotoharaetal.ll2002l) for TXS 0145+386 ER1 and ER2, as 
well as 4C 15.55 ER2. Deep J-band photometry with CISCO 
was also obtained for TXS 1812+412 ER2. All 6 OGs 
were observed in the i?-band using the 2048 x 2048 Tek- 
tronix (Tek2048) CCD on the UH 2.2m telescope, the Lo w- 
Resolution Imaging Spectrograph (LRIS: IOke et al.lfl995h on 
Keck I or K eck II, and/or the E chellette Spectrograph and 
Imager (ESI: ISheinis et aLll2002l) on Keck II. /-band imag- 
ing was also obtained with the Tek2048 CCD for the fields 
of TXS 1812+412 and TXS 0145+386. In addition, a deep 
z-band image of the field of TXS 0145+386 was obtained 
with LRIS and an RG850 filter. The combined RG850 filter 
with the LRIS CCD response is very similar to the standard 
SDSS z filter profile. Deep imaging was obtained for four 
of the sources (in the fields of B2 1018+34, TXS 121 1+334, 
4C 15.55, and TXS 1812+412) in both g- and R s filters on 
LRIS. The R s fil ter (also referred to as "sc ript R") is dis- 
cussed in detail by Steidel & Hamilton (1993); the R s magni- 
tude can be converted to a standard Cousins R using the R — I 
color index. TableQ]summarizes all of the observations. 

Flux cali bration was carried out using Landolt standard 
star fields dLando lt 1992) for the optic al imaging and UKIRT 
faint standards (Ha warden et al.ll2001l) for the near-IR imag- 
ing. The z-band imaging was an exception: because adequate 
z-band standards were not available at the time the observa- 
tio ns were taken, spectrophoto metric standards from the list 
oflMasseyi&Gronwall (1990) were observed instead. The 
magnitude zero point (Vega system) was determined by con- 
volving the Vega model SED with the normalized filter + 
CCD response function. 

The data were reduced using standard procedures in IRAF 
and photometry was performed on the resulting images us- 
ing apertures that were approximately 4" in diameter, the 
exact aperture size varying from source to source depend- 
ing on the extent of each individual galaxy. Aperture cor- 
rections were estimated from isolated stars in the field in or- 
der to match the 10" apertures that were used on th e stan- 
dard stars. We used values from Schleg eTet al.l (TT9 98) to cor- 
rect for Galactic extinction. Random errors were determined 
from the background variance and an additional 5% error was 
added in quadrature in order to conservatively account for ab- 
solute photometric calibration and systematic errors between 



the ground and space-based observations. All photometry is 
given in the Vega magnitude system. 

2.2. HSTACS andNICMOS imaging 

We also obtained HST imaging for all six sources. De- 
tails of the data reduction will be described in a future pa- 
per dMcGrath et al.l 120071) . Photometry was performed on 
the ACS F814W images using aperture radii that ranged be- 
tween l'/25 and l'/5. These values were chosen to include as 
much light from the galaxies as possible, while minimizing 
errors due to background noise. Preliminary modeling of the 
galaxy light profiles confirms that ^90% or more of the to- 
tal galaxy light falls within the measurement radii; any light 
falling outside the aperture is within the photometric errors. 
These aperture magnitudes were corrected to a 5" radius "infi- 
nite" aperture using correction factors estimated from isolated 
stars in the ACS field. The corrections were found to be in 
good agreement with the values given by Siri anni et al.l d2005l) 
based on PSF modeling. Finally, the aperture corrected instru- 
mental magnitudes were converted to Vega ma gnitudes us- 
ing the zero point for the F8 14W filter given by Sirian nTet al.1 

(2005) . 

NICMOS F160W photometry was performed in much the 
same way as that for the ACS images, however ap erture cor- 
rections were estimated from Tiny Tim (Rrist 2004) generated 
PSFs to determine what percentage of light from the PSF fell 
outside our measurement aperture. Galactic extinction cor- 
rections were applied to the HST data, and errors were deter- 
mined based on the rms background variation. A correction 
factor was applied to the rms errors for correl ated noise due to 
the drizzling process dFruchter & Hookl20 02). and a 5% abso- 
lute photometry error was added in quadrature to the random 
errors. Photometry from both ground and space-based obser- 
vations are given for the six sources in Table [2] The observed 
broad-band spectral energy distributions (SEDs) along with 
the best fit models are shown in Figure [2] and are discussed in 
detail in Sj5] 

2.3. Spectra of Five Old Galaxies 

Medium resolution spectra were obtained for five of the 
six old galaxies in our sample, with eith er LRIS on Keck 
I, or ESI or DEIMOS dFaberet al.l l2003h on Keck II. The 
LRIS spectrum of TXS 0145+386 ER1 and the ESI spec- 
trum o f 4C 15.55 ER2 are discussed in detail in lStockton et al.1 

(2006) . Spectra for B2 1018+34 ER2, TXS 1211+334 ER1, 
and TXS 1812+412 ER2 were obtained with DEIMOS on 
May 3-4, 2005. We used the gold-coated 1200 lines mm -1 
grating with a central wavelength of ~8000 A and a slit width 
of 1", which provided a spectral resolution of ~1.1 A. We 
used slits that were tilted by 5 degrees in order to improve 
sky subtraction. To minimize the effects of bad pixels and 
other chip defects, we designed the slit masks so that we could 
dither along the slit and maintain alignment. This was ac- 
complished through the use of small 1" alignment boxes at 5 
different dither positions 5" apart along the slit position an- 
gle and corresponding to the position of a faint star near the 
primary target. We used this faint star as a quick alignment 
check for each dither position. The data were reduced using 
the DEEP2 pipeline 4 and the resulting 1-d spectra for each in- 
dividual dither position were then co-added using the IRAF 

4 The spec2d analysis pipeline was developed at UC Berkeley with 
support from NSF grant AST-0071048 and can be obtained from 
http://astro.berkeley.edu/~cooper/deep/spec2d/ 
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TABLE 1 
Summary of Observations 



Field" 


Date 
(UT) 


Telescope 


Instrument 


Filter 


Exposure 
(s) 


TXS 0145+386 


1999 Oct 2 


Keck II 


LRIS 


Rc 


6300 




1999 Oct 12 


UH2.2 


Tek2048 


Ic 


20400 




1999 Oct 14 


UH2.2 


Tek2048 


Ic 


20400 




1999 Oct 2 


Keck II 


LRIS 


z' 


13500 




2005 Aug 16 


Subaru 


CISCO 


J 


960 




2005 Aug 16 


Subaru 


CISCO 


K' 


576 




2005 Jan 9 


HST 


ACS 


F814VK 


5288 


ER1 


2005 Jul 16 


HST 


NICMOS 


F160W 


2687 


ER2 


900S Till \(\ 

ZUUJ J 111 1 o 


HST 


1 N IV^IVIL/O 




2880 


ER1 


1999 Nov 3 


Keck I 


LRIS 


Spec 


16800 


ER1 


2001 Aug 24 


Keck I 


LRIS 


Spec 


10800 


B2 1018+34 


2005 Apr 9 


Keck I 


LRIS 


9 


5100 




2001 Apr 29 


UH2.2 


Tek2048 


Rc 


1800 




2005 Apr 9 


Keck I 


LRIS 


Rs 


5100 




2001 Feb 22 


IRTF 


NSFCam 


J 


1080 




2001 Feb 21 


IRTF 


NSFCam 


K' 


540 




2005 Jun 8 


HST 


ACS 


F814W 


5234 


ER2 


2005 Jun 8 


HST 


NICMOS 


F160W 


2688 


ER2 


9005 Mnv ^ A 


Kf*rV TT 




Spec 




TXS 1211+334 


2005 Apr 9 


Keck I 


LRIS 


9 


5440 




2002 May 18 


UH2.2 


Tek2048 


Rc 


5400 




2005 Apr 9 


Keck I 


LRIS 


Rs 


5100 




1998 Feb 13 


IRTF 


NSFCam 


J 


1080 




1998 Feb 14 


IRTF 


NSFCam 


J 


1080 




1998 Feb 12 


IRTF 


NSFCam 


K' 


540 




1998 Feb 14 


IRTF 


NSFCam 


K' 


540 




2005 Dec 4 


HST 


ACS 


FS1AW 


5234 


ER1 


2005 Jul 2 


HST 


NICMOS 


Fimw 


2688 


FJP 1 
1 IS. 1 


ZUUJ Ivldy J— 1 


J\.cLK 11 




Spec 


1 1 1 no 


4C 15.55 


2005 Apr 9 


Keck I 


LRIS 


9 


5780 




2002 Aug 7 


Keck II 


ESI 


R 


2520 




2005 Apr 9 


Keck I 


LRIS 


R s 


5700 




1999 Apr 8 


IRTF 


NSFCam 


J 


3780 




2002 May 31 


Subaru 


CISCO 


J 


960 




1999 Apr 8 


IRTF 


NSFCam 


K' 


2160 




2002 May 31 


Subaru 


CISCO 


K' 


1920 




2005 Feb 28 


HST 


ACS 


FS14W 


5000 


ER2 


2005 Mar 3 


HST 


NICMOS 


Fl&OW 


2688 


ER2 


2000 Jun 5 


Keck II 


ESI 


Spec 


18000 


TXS 1812+412 


2005 Apr 9 


Keck I 


LRIS 


9 


5440 




2002 May 18 


UH2.2 


Tek2048 


Rc 


1200 




2005 Apr 9 


Keck I 


LRIS 


Rs 


5100 




2001 Apr 29 


UH2.2 


Tek2048 


Ic 


3240 




2001 Apr 12 


IRTF 


NSFCam 


J 


1620 




2005 Aug 16 


Subaru 


CISCO 


J 


1920 




2001 Apr 12 


IRTF 


NSFCam 


K' 


1620 




2005 Jan 7 


HST 


ACS 


F814W 


5370 


ER2 


2005 Jun 8 


HST 


NICMOS 


F160W 


2688 


ER2 


2005 May 3-4 


Keck II 


DEIMOS 


Spec 


12560 



a The imaging fields were generally centered on the quasar, except for the HST NIC- 
MOS observations. Because of the small field of 19", these were centered on the OGs 
themselves, as indicated in column 1. For the HST ACS imaging, the quasars were 
placed near the center of the WFC1 CCD. The slit spectra were centered on the OGs. 



TABLE 2 
Photometry of Old Galaxies^ 



Galaxy 


g 


R 


I 


F814W 


z 


J 


F160W 


K' 


TXS 0145+386 ER1 




25.08±0.10 


23.28±0.()8 


23.26±0.09 


22.34±0.11 


20.54±0.08 


19.62±0.09 


18.51±0.06 


ER2 




25.52±0.14 


23.38±0.08 


23.21±0.08 


22.38±0.11 


20.64±0.08 


19.83±0.09 


18.75±0.08 


B2 1018+34 ER2 


25.81±0.11 


24.77±0.09 




23.00±0.08 




20.16±0.11 


19.59±0.10 


18.57±0.10 


TXS 1211+334 ER1 


27.14±0.22 


24.84±0.09 




23.02±0.08 




20.30±0.18 


19.26±0.07 


18.18±0.16 


4C 15.55 ER2 


> 26.96 (3<t) 


24.69±0.11 




22.80±0.06 




20.41±0.06 


19.42±0.08 


18.50±0.06 


TXS 1812+412 ER2 


26.74±0.25 


24.70±0.09 


23.04±0.13 


22.81±0.06 




20.47±0.08 


19.41±0.08 


18.36±0.06 


t Values are given in the Vega magnitude system and are corrected for Galactic extinction. 
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FIG. 2. — Broad-band photometry overplotted with solar metallicity instantaneous burst models from Chariot & Bruzual 12007). The best fit models both with 
(red) and without dust (blue) are shown. The values listed in the figure correspond to age in Gyr, magnitudes of rest-frame visual extinction, and P(x 2 ), the 
goodness of the fit. The model parameters are also given in Table[3] 
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task scombine. We cleaned and smoothed the spectra to a 
resolution of ~8.5 A in order to more easily detect spectral 
features. We also obtained slitless spectroscopy of Feige 34 
at the position of each primary target in order to flux cali- 
brate the spectra, while t he standard Mauna K ea models of 
atmospheric extinction (Krisciunas et al1 ll987l) were used to 
correct for small (±0.3) airmass differences. Finally, the flux 
calibrated spectra were scaled by small amounts to account 
for slit losses and to match the R- and /-band photometry. 
The final spectra are shown in Figure [3] al ong with the best fit 
prelim inary stellar population models of ICharlot & Bruzuall 
(2007) 5 . 

For all five sources, we were able to determine redshifts 
from either absorption or emission lines in the spectra. In 
four out of five cases, these redshifts are close to that of the 
radio QSO (Az < 0.01). For TXS 1812+412 ER2, the red- 
shift is significantly different from that of the radio source. 
The galaxy therefore cannot be associated with the quasar en- 
vironment and is, rather, a chance projection. Spectroscopic 
redshifts are given along with the QSO redshifts in Table [3] 

3. GALAXY SPECTRAL ENERGY DISTRIBUTIONS 
AND AGES 

For our analysis, we have used ICharlot & Bruzuall d2007l 
CB07) instantaneous burst models with a IChabrieri (120031) 
IMF, solar metalli city, and the empirical Calzetti law 
dCalzetti et al.l 120 00) for dust extinction. Broad-band photo- 
metric fits to the mo del SEDs were com puted with a version 
of the Hyperz code {Bolzonella et al. 2000) that was modified 
in order to allow finer age steps. The fits for all six of the OGs 
are shown in Figure [2] For all cases except TXS 0145+386 
ER2, the redshift has been fixed to the spectroscopically de- 
termined value so that age and reddening are the only free 
parameters in the fit. For TXS 0145+386 ER2, the redshift 
range was allowed to vary around the values for the QSO and 
that of ER1, and the best-fit value was found to be z = 1.459. 

It is important to note that there is quite a bit of uncer- 
tainty in the stellar population models themselve s, as well as 
the tr eatment of dust for high redshift sources dPierini et al.l 
|2005l) . In addition, thermally pulsating asymptotic giant 
branch (TP-AGB) stars are known to contribute a significant 
amount of flux in the inf rared for stellar populations of in- 
termediate age (~1 Gyr; Maraston 2005). The amount of 
this contribution is not well known; however, the net effect 
of including TP-AGB stars is that the population synthesis 
models become redder overall, allowing younger, less dusty 
fits to observed galaxy SEDs. The CB07 models include a 
larger contribution fro m TP-AGB stars than any other models 
to date, including the Maraston] (120051 M05) mod els, result- 
ing in younger age determinations (Bruzua fl2007h . In prac- 
tice, however (as discussed in fQ, we find that by and large 
when dust is included, the various models all tend to converge 
towards similar ages. We therefore believe that by using these 
improved models, our derived galaxy ages and formation red- 
shifts should be fairly reliable estimates. 

Another source of uncertainty c omes from the well-known 
age-metallicity degeneracy (e.g., IWorthevl fl994). On aver- 
age, a factor of 2.5 increase in the metallicity of the model 
stellar population results in a best-fit age that is 30-50% 

5 These models are preliminary in the sense that they do not yet include the 
updated Padova 2007 ste llar evolutionary tracks, although they do include the 
Marigo & Girardi 12007) prescription for thermally-pulsating asymptotic gi- 
ant branch evolution of low- and intermediate-mass stars (S. Chariot, private 
communication). 



younger. An important caveat to this simplistic view of the 
age-metallicity degeneracy is that increasing metallicity can 
only be balanced by younger ages to a certain point. Ages 
much younger than ~1 Gyr do not provide good fits to the 
data regardless of metallicity as they simply cannot repro- 
duce the sharp inflection at (rest-frame) 4000 A required by 
the observations. Furthermore, the ages are determined to a 
large extent by both the 4000 A break and the spectral slope 
in the rest-frame near-UV. For realistic chemical evolution 
models, this spectral region is likely to be dominated by sub- 
solar-metallicit y stars even if the average stellar metallicity 
is super-sol ar (Moll enFritze-v. Alvensleben. & Frickdll997t 
iBicker et al.l2004l) . Briefly, this result follows from the spre ad 
in metallicities that will always be present in the stellar pop- 
ulations in a galaxy and the lower line blanketing in the near- 
UV for lower- metallicity stars. 

From Figure [2] and Table [3] it is clear that the majority of 
these z ~ 1.5 galaxies are dominated by old stellar popu- 
lations, implying that star formation ceased at even higher 
redshifts. SED fits that do not include dust yield older ages; 
however, in many instances our current photometry is insuf- 
ficient to distinguish between the best fit dust-reddened and 
no-dust models. The difference between these models often 
amounts to only a few percent change in the overall quality 
of the fit, such that the no-dust model cannot be ruled out 
with any certainty. It is also important to note that, while the 
fits may be marginally better including dust, to some extent 
this is expected solely because there are more free param- 
eters in the fit. Longer wa velength data from Spitzer (e.g., 
IStockton & McGrafhli2007l) is needed to clarify how much 
dust is present in these galaxies and to better constrain their 
ages. 

For the most part, however, the SEDs paint a consistent pic- 
ture. With the exception of B2 1018+34 ER2, the galaxies are 
dominated by stellar ages that are ~ 1 Gyr or greater. The most 
extreme case is TXS 0145+386 ER2, which is dominated by 
a ~2.5 Gyr stellar population, with little or no dust. 4C 15.55 
ER2 also seems to be entirely free of dust, with an age of 1.8 
Gyr. Assuming purely passive evolution following a single 
massive burst of star formation, the implied formation red- 
shifts for these two galaxies are Zf — 3.3 and Zf = 2.4, 
respectively. More complex and more realistic star formation 
histories, such as exponentially decaying or truncated star for- 
mation, would imply even earlier average formation times. 

On the younger edge of the age distribution, it is possi- 
ble to fit the photometry for both TXS 0145+386 ER1 and 
B2 1018+34 ER2 with models where the break between I- 
and J-bands is due to the heavily reddened Balmer break of an 
intermediate age stellar population instead of a strong 4000 A 
break from older stellar populations. In fact, B2 1018+34 ER2 
was inadvertently allowed a less stringent criterion for J—K 1 
during the initial selection, and we can see clearly that this 
galaxy does not fit the same mold as the other galaxies. A 
significant amount of dust (with a younger underlying stel- 
lar population) appears to be almost required in order to fit 
the bluest (g-band) data point, and the J-band point is likely 
contaminated by an emission line due to star formation or an 
active nucleus. H/3 falls within the J-band at z = 1.4; how- 
ever, the lack of any apparent excess in the NICMOS F160W 
photometry, where we would expect to see Ha, rules out H/3 
emission sufficiently strong to produce the observed excess 
in the J band. This excess could, however, be due to strong 
[O III] emission. The spectrum of this galaxy (Figure |3]l does 
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FIG. 3. — The rest-frame UV spectra plotted in red for five of the six galaxies in our sample. Overplotted in black is the best-fit dust-reddened Chariot & Bruzual 
12007) solar metallicity, instantaneous burst model determined from the photometric SED fitting. The resolution of the model at A < 3300 A is ~ 20 A, while at 
A > 3300 A, it has been smoothed to ~ 8.5 A in order to match the observed data. On the left, TXS 0145+386 ER1 and B2 1018+34 ER2 show [O II] emission, 
which could be due to either recent star formation or active nuclei. On the right we show the three oldest galaxies. The Mg II 2800 A feature is clearly visible 
in TXS 121 1+334 ER1 and 4C 15.55 ER2, while Ca II H and K are visible in TXS 1812+412 ER2. Residual A-band atmospheric absorption in the spectrum of 
TXS 121 1+334 ER1 is marked for reference. The ~300 A-wide region of the spectrum of TXS 1812+412 ER2 around rest-frame 3550 A is an artifact. Redshifts 
determined from absorption and/ or emission features in each galaxy's spectrum are given in Table|5] 



show [O II] emission, and we disc uss t his, as well as possible 
[O III] emission in more detail in 1 1 

A dust reddened model seems favored for TXS 0145+386 
ER1, as well, although the underlying stellar population ap- 
pears to be older than that of B2 1018+34 ER2. A deep blue- 
band photometry point or deep mid-infrared data would help 
clarify the amount of dust present, and truly pin down the age 
of the galaxy. Over the observed wavelength range, there is 
very little difference between the dust-reddened and no-dust 
models. Given the dust-free model, the implied formation 
redshift for this source is z — 3.7. If we assume, however, 
that the dust-reddened model is correct (as it exhibits a higher 
probability than the no-dust model) the implied formation red- 
shift for this source is a modest zt = 1.9. The spectrum 



(Figure [3} shows that [O II] is present in this galaxy, as well, 
lending credence to the idea that a fair amount of dust could 
be present. This [O II] emission was discussed in detail in 
Stock ton" et al.l (120061) . and is either due to a weak AGN or 
a starburst-driven wind since the linewidth implies a veloc- 
ity that is too great to be due to ordered rotation or internal 
velocity dispersion alone. 

The remaining two sources, TXS 1211+334 ER1 and 
TXS 1812+412 ER2, are well-fit by intermediate amounts of 
dust (0.23 and 0.51 magnitudes visual absorption, respec- 
tively) and underlying old stellar populations. The difference 
in x 2 between the dust-reddened and no dust models, while 
not overwhelmingly large, does favor the presence of at least 
some dust in these galaxies. In fact, given the presence of 
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TABLE 3 

Redshifts and Model Parameters 



Galaxy 


z QSO 


Zspec 


Model 


Age— no dust 


Age— with dust 


Ay 


Mass 










(Gyr) 


(Gyr) 


(mag) 


(10 n M Q ) 


TXS 01 45+386 FR 1 


1.442 


1.4533 


CB07 


2.75 


0.90 


0.79 


1.74 








BC03 


3.00 


1.28 


0.57 


2.00 


ER2 


1.442 


1.459* 


CB07 


2.50 


2.40 


0.02 


2.02 








BC03 


2.75 


2.50 


0.08 


2.27 


B2 1018+34 ER2 


1.4058 


1.4057 


CB07 


1.14 


0.64 


0.94 


1.42 








BC03 


1.14 


0.57 


1.17 


1.68 


TXS 1211+334 ER1 


1.5960 


1.598 


CB07 


1.68 


1.43 


0.23 


2.90 








BC03 


1.80 


1.43 


0.29 


3.15 


4C 15.55 ER2 


1.406 


1.412 


CB07 


1.80 


1.80 


0.02 


1.75 








BC03 


2.10 


1.90 


0.09 


2.06 


TXS 1812+412 ER2 


1.564 


1.290 


CB07 


3.00 


1.43 


0.51 


1.56 








BC03 


3.50 


1.28 


0.73 


1.79 



* The best-fit photometric redshift is given as no spectroscopy was obtained for this source. 



[O II] in the previous two galaxies (B2 1018+34 ER2 and 
TXS 0145+386 ER1) and the ability of the population syn- 
thesis models to reproduce a moderate amount of dust from 
the photometric fitting independently of the spectra, gives us 
confidence that the dust-reddened models provide reasonably 
good representations of the galaxy SEDs. It is important to 
note, however, that models which imply both large inferred 
ages and high global extinction values would be difficult to 
explain from a physical standpoint given the lack of evidence 
for any significant amounts of recent star formation. In these 
cases, further observations combined with better modeling 
would help clarify the situation. With dust, the implied for- 
mation redshifts for TXS 1211+334 ER1 and TXS 1812+412 
ER2 are Zf — 2.5 and Zf = 1.9, respectively. Without dust, 
these values shift to Zf = 2.7 and Zf — 3.4. 

4. STELLAR POPULATIONS FROM REST-FRAME UV 
SPECTROSCOPY 

In the following, we discuss the spectra of the three sources 
we observed with DEIMOS on Keck II. For a detailed dis- 
cussion of the LRIS spectrum of TXS 0145+386 and the ESI 
spectrum of 4C 15.55 ER2, refer to lStockton eTaT] d2006l) . We 
have re-fit the spectra for TXS 0145+386 ER1 and 4C 15.55 
ER2 using the CB07 models determined from the photometry 
in order to be consistent with the remaining three sources, and 
these are shown together in Figure [3] 

4.1. [O II] emission in B2 1018+34 ER2 

The spectrum of B2 1018+34 ER2 exhibits [O II] emis- 
sion (3727 A, Figure [3]), pointing to the possibility of recent 
star formation in this galaxy. The linewidth for the [O II] 
emission is a ~ 200 km s _1 , which is consistent with ro- 
tational broadening or internal velocity dispersion. If we as- 
sume the [O II] emission is due to recent star formation, the 
equivalent width of the [O II] line implies a star formation 
rate (SFR) of ~ 3.5Mm yr" 1 , us ing the relationship given in 
I Kewlev. Geller. & Jansen] d2004l eq. 19), which is calibrated 
to agree with Ha SFR indicators and is corrected for extinc- 
tion by assuming Ay ~ 1. Given the best fit SED to the 
broadband photometry (FigureO, it is likely that a significant 
amount of dust, up to 1 .0 magnitudes rest-frame visual extinc- 
tion, is present in this galaxy, indicating that this is a reason- 
able estimate of the star formation r ate. This SFR is simil ar to 
that found for TXS 0145+386 ER1 (Stockto nltaill200^) and 
to that of dusty, star-forming EROs in the lCimatti et al.l(l2002l) 



sample. On the other hand, if the observed excess in J band is 
due to [O III] emission, the implied [O lll]/[0 II] ratio would 
be —24, which is too great to be explained by anything other 
than an AGN. Such strong [O III] emission would easily be 
detectible with J-band spectroscopy, which might also have a 
good chance of detecting broad H/3, if present. 

Regardless of whether the [O II] emission is due to star for- 
mation or an AGN, it is clear that the underlying stellar pop- 
ulation is still relatively old. A CN absorption feature is vis- 
ible at 3831 A as well as a hint of the Mg II line at 2800 A; 
however, the latter is very near the edge of the spectral sensi- 
tivity. Assuming a single stellar population, the implied for- 
mation redshift for this source (zf ~ 1.7) is easily within the 
allowable range of most if not all current SAMs of galaxy for- 
mation. If there is ongoing star formation, it is a relatively 
small effect by mass, while if there is an AGN present, then 
it could be the catalyst that quenched the star formation. The 
timescale for the latter case app ears to be in line wit h the AGN 
feedback scenario proposed bv lCroton et al.l d2006l) . 

4.2. Spectrum of TXS 1211 +334 ER1 

The DEIMOS spectrum of TXS 1211+334 ER1 (FigureE) 
shows a very clear absorption trough with a shape character- 
istic of the Mg II 2800 A line. The strength of this absorption 
line, as well as the slope of the continuum are well fit by a 
1.43 Gyr CB07 solar metallicity instantaneous burst model 
with 0.23 magnitudes of extinction, which was determined 
through the photometric SED fitting. A least squares fit to 
the spectrum alone yields an age of 1.9 Gyr, although the chi- 
squared has a shallow minimum from 1 .5 to 2.0 Gyr. The red- 
shift derived from the Mg II line is z = 1.598, which when 
combined with the age from the photometry implies that by 
z = 2.5 the galaxy had completely ceased star formation and 
had accumulated nearly 3 x lO n -M0 (Table|3]l. 

4.3. Spectrum of TXS 1812+412 ER2 

The spectrum for TXS 1812+412 ER2 shows a broad ab- 
sorption feature near 8100A(rest-fr ame -3550 A). This fea- 
ture cannot be attributed to atmospheric absorption, both be- 
cause of the width of the observed feature (—300 A observed 
frame) and because a secondary source on the slit roughly 10" 
away does not show a similar drop in intensity. The feature is 
also visible in multiple dither positions. However, the broad- 
band photometry (both /-band and F814W observations) do 
not show any signatures of this absorption. We therefore are 
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confident that this feature is an instrumental or data-reduction 
artifact, and we refrain from using this portion of the spectrum 
in our analysis. 

The best-fit CB07 model from the photometry of a 1.43 
Gyr stellar population with 0.51 magnitudes of extinction is 
shown overplotted on the observed spectrum (FigurefJ]). Note 
that while many features, such as Ca II H and K, and CN at 
3831 A, are fairly well reproduced by this model, the depth of 
the lines and the underlying slope of the continuum are not. 
This is suggestive of either an older age or a higher metal- 
licity stellar population. Recall that we restricted our photo- 
metric fits to solar metallicity models. A 2.5 x Z Q model 
with an age of 0.72 Gyr plus 0.76 magnitudes of visual ex- 
tinction fit to the broad-band photometry provides a slightly 
better fit to the overall spectrum, as does the dust-free 3.0 Gyr 
solar metallicity model. Fitting the spectrum alone, without 
including longer wavelength photometry, results in a best-fit 
2.5 x Zq model of 1 .4 Gyr, while a Zq model yields an age of 
~3.5 Gyr. Given that there are no signatures of young stellar 
populations or post-starburst features in the spectrum, older 
ages are slightly favored, however an age of 0.72 Gyr is not 
necessarily precluded by this, especially if moderate amounts 
of dust are present as inferred from the best fit photometric 
model. Once again, longer wavelength data is needed in order 
to better discriminate between models. 

5. DISCUSSION 

Previous st udies of passive ga l axies a t z > 1 have largely 
relied on the iBruzual & Chariot) (I2003L BC03) stellar popu- 
lation synthesis models, which include a much lower con- 
tribution from TP-AGB stars than do th e CB07 models, for 
estimates of galaxy age and mass (e.g., ICimatti et alj 



esti mates ot galaxy age and mass (e.g., ICimatti et a. 
McC arthy et al.1 12001 lYan et al l 120041: iDaddi et al 



2004 



Saracco et all 120051: iKriek et al.l 120061; IStockton et all 



2005 



2004 



2006). Others have relied on the M05 models (Pieri ni et al.l 
20051: iMaraston et al .1120061) . which are generally more simi- 
lar to the CB07 models with respect to the TP-AGB evolution- 
ary phase. In order to compare our results with these studies, 
we have calculated the best-fit BC03 model SED fits to our 
photometry, as well as the best-fit BC03 and M05 models to 
the galaxy spectra. The BC03 model fits to the photometry 
are presented in Table [3] along with the best-fit CB07 mod- 
els. When reddening is allowed to be a free parameter, in 
all but one case (TXS 0145+386 ER1) the BC03 models re- 
quire more dust to fit the photometry. If we hold E(B-V) = 
0, the best-fit BC03 models are older than the CB07 no-dust 
models for five of the six galaxies, but not dramatically so. 
The average increase in age is ~10% for the BC03 models. 
When we include dust, the difference in age between the two 
models is negligible, and does not show any preferred trend. 
This implies that TP-AGB stars do not greatly influence the 
shape of the SEDs for these galaxies, perhaps because the 
age is old enough that TP-AGB stars have all evolved to the 
point where their contribution to the total bolometric lumi- 
nosity is dwarfed by lower mass, main sequence stars. In 
addition, the rest-frame near-IR, where the contribution from 
TP-AGB stars is greatest, is not explored by our observations. 
The clearest difference between the BC03 and CB07 models 
is that our inferred galaxy masses are slightly smaller using 
the CB07 models (^10%), and in general we find that the \ 2 
fits are better with the CB07 models than with the BC03 mod- 
els. 

When we constrain our model fitting to the spectra alone, 
we get results that are typically within a few hundred Myr of 



the values determined using similar models with the broad- 
band photometry. The only exception to this is that for 
4C 15.55 ER2, the best-fit CB07 age from the spectrum is 1.1 
Gyr, as opposed to an age of 1.8 Gyr from the photometry. 
The 1 . 1 Gyr model, however, provides an extremely poor fit 
to the near-IR photometry, with a P(x 2 ) close to 0%, while 
the chi squared fit to the spectrum has a broad minimum from 
1.0 to 1.8 Gyr. The older age from the photometry is there- 
fore preferred, as it is the most consistent with both datasets. 
When fitting the spectra, the difference in x 2 between dust 
and no-dust models is often very small, as would be expected 
given the small wavelength range of the fit; the presence of 
dust effectively just changes the scaling required to fit the 
data. This is exemplified in the fit for TXS 0145+386 ER1, 
where the difference in \ 2 between the no-dust (2.75 Gyr) 
and dust-reddened (1.1 Gyr) models amounts to only a 7% 
difference in the overall quality of the fit. We find that the 
M05 models agree well with the CB07 models, however the 
average best-fit age is slightly older when using the M05 mod- 
els, and the agreement with the near-IR photometry is gener- 
ally worse. Furthermore, when we allow metallicity to vary, 
the best-fit models favor solar-metalli city in all but one case 
(TXS 1812+412 ER2, as discussed in gjj. 

The seven passive galaxies found i n the Hubble Ultra 
Deep Field bvlDaddi et alj (120001 [2005) and re-examined by 
Maras ton et alJ(l2006l) . have similar formation redshifts to five 
of our six z = 1.5 galaxies and therefore provide a fair com- 
parison to the galaxie s in the present study. Passive galax- 
ies in the IDaddi et aH sample are generally a factor of two 
less massive than ours and a few exhibit a brighter blue tail 
or bump in the SED which yield ages as young as 0.3 Gyr. 
Several of their objects at z < 2, however, have properties 
similar to the ones in this study, which expands on the num- 
ber of known high reds hift evolved galaxi es. At even higher 
redshifts (2 < z < 2.7). lKriek et al.l d2006l) discuss nine mas- 
sive galaxies whose stellar continua and Balmer breaks are 
detected in near-IR spectroscopy, but show no trace of emis- 
sion lines from young star formation. These galaxies are fit 
by stellar populations with ~1.0 magnitudes visual extinction 
that are on average ^500 Myr old, which imply formation 
redshifts of Zf ~3. The galaxies in the present study could 
very well be the passively evolved descendants of such a pop- 
ulation, wit h the proviso that much of the dust present in the 
IKriek et ail sample would need to be expelled rather than con- 
sumed in star formation between 2.5 < z < 1.5 in order to 
match both the lower dust content and the continued passive 
evolution seen in our z = 1.5 galaxy sample. As a caution- 
ary note, however, the ages for these 2 < z < 2.7 galaxies 
fall within the range where the contribution from TP-AGB 
stars may be much more significant than for our ~2 Gyr-old 
z = 1.5 galaxies, and therefore care must be taken in com- 
paring the two populations. Use of models with improved 
treatments of TP-AGB stars such as M05 or CB07 could yield 
lower ages for these sources. Alternatively, the M05 or CB07 
models could require less dust, thereby bringing these obser- 
vations into even better agreement with the idea that they are 
progenitors to galaxies at z = 1.5 both in the present study 
and in the Dad di et all (120051) sample. Finally, galaxies such 
as the evolved disk at z = 2.5 found bv lStockton et al.l (|2004|) 
appear to be altogether a much more extreme case of early 
and rapid formation of the most massive galaxies, with an im- 
plied formation < 1 Gyr after the Big Bang. It is unclear how 
common objects such as these are, and it will be important to 
continue searching for the oldest stellar populations at higher 
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redshifts in order to further constrain galaxy formation mod- 
els. 

6. CONCLUSIONS 

We have studied a sample of six evolved galaxies at z ~ 
1.5 with broadband photometry covering the rest-frame UV 
to visible wavelengths, as well as medium resolution spec- 
troscopy from Keck covering the rest-frame UV portion of 
the spectrum. Us ing the most recent stellar population syn- 
thesis models of ICharlot & Bruzuall d2007j) which include 
an improved treatment of TP-AGB stars, we find ages for 
these galaxies that imply formation redshifts roughly between 
2 < z < 4. Simple SAMs are unable to reproduce mas- 
sive galaxies at such high redshifts, however, models which 
employ feedback from AGN may not be as inconsistent with 
these observations. Galaxies such as TXS 0145+386 ER1 and 
TXS 1812+412 ER2, under the assumption of low dust con- 
tent, test the limit of galaxy formation theories by requiring 
massive galaxies to be present at z > 3.4. If dustier models 
are preferred for these galaxies, however, the formation red- 
shifts are more in line with values expected in the AGN feed- 
back scenario. Clearly, longer wavelength data from Spitzer 
is needed to pin down the dust content and age of these galax- 
ies. Irrespective of dust, the best-fit spectral synthesis model 
for TXS 0145+386 ER2 implies a 2 x 10 n M Q galaxy that 
completed star formation by at least z = 3.3. This may be 
difficult to explain even with AGN feedback. 

[O II] is observed in two of the sources (TXS 0145+386 
ER1 and B2 1018+34 ER2), which may either be an indica- 



tion of lingering star formation or the remains of weak AGN. 
If the latter scenario is correct, it may give an indication over 
what timescale AGN activity is observable after star forma- 
tion has ceased, which may prove useful for testing the AGN 
feedback scenario. 

Overall, these observations imply that massive galaxies did 
in fact exist at very early times in the universe. Even given 
the uncertainties in the population synthesis modeling, and 
the trend toward lower ages for models with improved treat- 
ments of TP-AGB stars, we are faced with several examples 
of massive galaxies (> 10 11 Mq) that appear to have been 
passively evolving since z > 2. We will exp lore how these 
galax ies may have formed in a future paper (M cGrath et alJ 
120071) through a high-resolution morphological study. 
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